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In this work we show that Rashba-based topological superconductor nanowires, where the spin- 
orbit coupling may change its sign, support three topological phases protected by chiral symmetry. 
When a superconducting phase gradient is applied over the interface of the two nontrivial phases, the 
Andreev spectrum is qualitatively phase shifted by tt compared to usual Majorana weak links. The 
topological 7r-j unction has the striking property of exhibiting maximum supercurrent in the vicinity 
of vanishing phase difference. Qualitative features of the junction are robust against disorder and 
magnetic fields violating chiral symmetry. The studied system could be realized by local gating of 
the wire or by an appropriate stacking of permanent magnets in synthetic Rashba systems. 
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Introduction- Recent experimental evidence is consis- 
tent with the possibility that superconducting nanowires 
with a strong spin-orbit coupling support topological su- 
perconductivity and Majorana bound states (MBS) [IH3]. 
Though the existence of MBS is not conclusively settled 
yet, experimental efforts have increased the interest to- 
wards the topic significantly. Besides being interesting 
due to their exotic physical properties, MBS could find 
applications in quantum information processing. [4-8J. 

Theoretical developments in topological insulator 
physics [9, 10J and two-dimensional semiconductor het- 
erostructures [11, 12J paved the way for the proposal of 
topological superconductivity in nanowires with a strong 
Rashba spin-orbit coupling p~3j Hi] , The topological 
phases of nanowire systems are commonly classified as- 
suming only particle-hole symmetry leading to a Z<i in- 
variant. However, in the special case when the applied 
magnetic field is perpendicular to the effective spin-orbit 
field and the superconducting order parameter is uni- 
form, the system possesses additional chiral symmetry. 
Then the system may have an arbitrary number of dis- 
tinct topological phases characterized by a Z-valued in- 
variant fT5UT7] , Since unbroken chiral symmetry relies on 
fine tuning of the magnetic field direction, classification 
of the nanowire system is usually regarded as Z2 type. 

In this work we show that the concept of chiral symme- 
try in nanowires is not merely a mathematical curiosity 
and can have distinct physical consequences. We con- 
sider a nanowire system where in some segments of the 
wire the Rashba coupling constant may have opposite 
signs and show that this system possesses three topo- 
logical phases protected by chiral symmetry. The three 
phases are distinguished by an invariant taking values 
v = ±1 and v = 0. The phases v = zbl correspond to 
the nontrivial phase in the Z2 classification while v — 
phase coincides with the trivial phase. The interface be- 
tween ±1 and is similar to the boundary between trivial 
and nontrivial phases in Z2 classification but a boundary 
between 1 and —1 phases contains two degenerate MBS 
protected by chiral symmetry. As discussed below, the 



sign change of the Rashba coupling could be realized, 
for example, by local gating of the wire. Alternatively, 
in synthetic Rashba systems based on spatially varying 
magnetization, an appropriate stacking of magnets in the 
proximity of the wire could implement the sign reversal. 

When a superconducting phase difference is applied 
over the phase boundary where the Rashba coupling 
changes sign, chiral symmetry is violated and the degen- 
eracy of the two MBS in the junction is lifted. The two 
MBS are degenerate when the phase difference is 2n7r, 
in contrast to usual Majorana weak links, spectrum of 
which has a pair of degenerate MBS when the phase dif- 
ference is (2n + 1)tt with an integer n. Thus, the sign 
change of the Rashba coupling qualitatively amounts to 
a phase shift of tt of the Andreev spectrum. Due to spe- 
cial properties of a topological Josephson junction, the 
spectrum of a topological 7r-j unction supports maximum 
supercurrent around vanishing phase difference. 

We begin by introducing the nanowire model and clas- 
sifying phases according to chiral symmetry. The MBS 
wavefunctions on the —1/1 phase boundary are solved 
analytically in special cases. The Andreev spectrum of a 
topological 7r-j unction is studied numerically by solving a 
tight-binding model of a finite nanowire, including impu- 
rity potential and magnetic fields violating chiral symme- 
try. Qualitative features of the junction are insensitive 
to the detailed manner how the Rashba coupling changes 
sign. 

Model and phase diagram- We study a wire with a 
Rashba spin-orbit coupling in the vicinity of a supercon- 
ductor in magnetic fields. The system is modeled by the 
Bogoliubov-de Gennes Hamiltonian p~3| H3], 

H(k) = (e k + aka x ) r z + B • a + A • r, (1) 

where e k = %^ ^ and cr^ and are Pauli matri- 
ces operating in the spin and the Nambu space, respec- 
tively. The Hamiltonian is expressed in the basis 
\I/ = — ipl) T and the coordinates have been 
chosen so that the wire is parallel to the y direction and 
the effective spin-orbit field coincides with the x direc- 
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Figure 1: a): Phase diagram of the nanowire system in 
(B 2 — A 2 — /i 2 , a)-plane in the presence of chiral 
symmetry, b): Schematic representation of the 
topological 7T-j unction formed between wire segments 
where the Rashba coupling changes sign. The interface 
contains two Majorana bound states that become 
(nearly) degenerate at the vanishing superconducting 
phase difference. Inset: Magnetic realisation of an 
effective sign-changing Rashba coupling. 



tion. The first term represents kinetic energy of elec- 
trons and holes including the Rashba coupling charac- 
terized by a coupling constant a. The second term is 
the Zeeman coupling B = (B x1 B yi B z ) due to the ap- 
plied magnetic field and the last term arises from the 
proximity-induced superconducting pairing with A = 
A(cos((p)r x , sm((p)r y , 0). 

In the special case when the magnetic field is per- 
pendicular to the spin-orbit field and the superconduct- 
ing phase is uniform in the system, one can choose 
B = (0,0, B z ) and A = (A, 0,0). Then the Hamilto- 
nian anticommutes {H, C} = with C = cr x T y , thus 
possessing chiral symmetry. As a consequence, if H 
has an eigenstate \Ef) with energy E^ then there ex- 
ists another eigenstate C\Ef) belonging to energy — Ei. 
Zero modes, when present, can be chosen to coincide 
with the eigenstates of the generator of chiral symmetry, 
C\^) = ±1^^). In odd spatial dimensions systems with 
chiral symmetry posses topological phases characterized 
by a Z- valued topological invariant [T5HT7] . There exists 
different equivalent formulations of the invariant, here we 
consider a form appropriate for Id systems 



u[H] = - 1 - [ dkTr ICH^dkH] 
Aiti J L J 



(2) 



which was recently employed in Ref. [18] to analyse sim- 
ilar insulating systems. The phase diagram of the sys- 
tem in the (B 2 — A 2 — /i 2 , a) plane is depicted in Fig. 1 
a), revealing three distinct phases v = ±1 and v — 0. 
The v — phase, corresponding to the parameter regime 
B 2 < A 2 + /i 2 , can be identified with the trivial phase 
in the Z 2 classification. The v = ±1 phases differ by the 
sign of the Rashba coupling and map to the nontrivial 
phase in the Z 2 scheme which is insensitive to the sign 
of a. Thus the phase boundaries between v = ±1 and 
v = phases behave exactly as the Z 2 interface between 
the trivial and the nontrivial phase supporting a single 



MBS [14J. However, the boundary between v = — 1 and 
v = 1 is expected to be qualitatively different and, in 
fact, it supports two MBS. These MBS are degenerate 
at zero energy and cannot be fused to form finite energy 
states as long as chiral symmetry remains intact. 

The existence of two MBS can be demonstrated by 
direct analytical solution of the zero modes on the —1/1- 
phase boundary. To simplify the analysis, we assume 
that the Rashba coupling a(x) is negative constant for 
x < and positive constant for x > and that the 
parameters satisfy /i, A < B <C Er = ma 2 /2h. Then 
we can linearize ([TJ in k and solve the zero mode wave 
functions from [—ia(x)a x r z d x — jir z + Ba z + Ar x ] ip = 
|14j . Matching the solutions for different signs of a at 
x = we recover two normalizable solutions 
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where e %e = ^ — i y ^ B - and Ai/ 2 are real constants 
(here B 2 > A 2 + /i 2 ). The vectors \x±) and \y±) 
correspond to the eigenstates of operators a x and r y 
and the upper (lower) signs on the right-hand side of 
Eq. ([3| correspond to ipi (^2)- The zero modes satisfy 
C'0i/2 = i'01/2 and can be obtained analogously for the 
opposite boundary where v — 1 changes to v — — 1. Ex- 
pressing the wavefunctions in a standard four component 
form, the corresponding field operators can be written as 
7i/2 = ^1/2 " ^5 which satisfy 71/2 = 7i/2- Although the 
analytical MBS solutions were obtained in a special pa- 
rameter regime, the number of solutions cannot change 
as long as the bulk gaps are finite and chiral symmetry 
remains unbroken. 

Topological tt -junction: general discussion- Different 
aspects of the Josephson effect in topological supercon- 
ductors are under active research fT9H26j . The 47r pe- 
riodic Josephson effect between two topological p-wave 
segments coupled through a weak link was already dis- 
cussed in the pioneering work by Kitaev |5J . Typically the 
Andreev spectrum of a topological nanowire weak link 
contains two MBS states that become degenerate when 
the superconducting phase difference cp = (2n + l)7r is ap- 
plied over the junction. This is qualitatively in striking 
contrast to the above discussed behaviour the v — ±1 
interface which supports degenerate MBS at ip = n2ir. 
This suggests that spectrum of the v = ±1 junction is 
qualitatively phase shifted by tt. 

As discussed by a number of authors, Hamiltonian 
can be expressed in the basis ip± of eigenfunctions of the 
system in the absence of superconductivity (see for exam- 
ple Refs. [27 ; 28J). At low filling only the lower band ip- 
is populated and the effective gap function projected to 
this band is given by A p = -y=^==A. This expression 
suggests that changing the sign of the Rashba coupling a 
is effectively equivalent to imposing a rigid phase shift of 
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tt to the order parameter. Therefore we expect that the 
phase boundary between the v = 1 and v = — 1 phases 
effectively behaves as a 7r-j unction when a finite phase 
difference is applied. The topological 7r-junction has an 
interesting difference compared to the usual tt-j unctions. 
In the usual case the tt phase shift replaces a local min- 
imum of the original spectrum at tp = (mod 2tt) by 
a local maximum. Since the current through a short 
junction is given by I((p) = ^d^E^ip), where E(cp) is 
the energy of the populated Andreev level, supercurrent 
vanishes at ip = for a topologically trivial tt-j unction. 
In the topological 7r-j unction the tt phase shift moves the 
zero-energy crossing from (p = tt to ip = (mod 2tt). The 
vicinity of the crossing point is not an extremum of the 
spectrum, therefore giving rise to a spontaneous super- 
current 1(0) 7^ in infinite systems. In a finite system 
the picture is slightly modified since the crossing at (p = 
becomes avoided due to coupling to the MBS at the end 
of the wire (depicted in Fig. 1 b) ) and supercurrent van- 
ishes 1(0) = 0. Avoided crossing has the same origin as 
in usual topological Josephson junctions [20, 21J. How- 
ever, the avoided crossing becomes exponentially narrow 
when the wire is longer than the localisation length of 
MBS wavefunctions, thus yielding maximum supercur- 
rent in the immediate vicinity of cp = 0. 

Tight-binding calculation of tt -junction spectrum- 
Above we gave general arguments why we expect v = ±1 
interface to behave as a 7r-j unction. Here we study spe- 
cific properties of the topological 7r-j unction by numeri- 
cally solving the spectrum of a tight-binding model of the 
setup in Fig.[l]b). The Hamiltonian of a lattice model is 
given by 

H = -J2( t + M + ^(-/itfKtfi + V^lr^,) 

+ ^ Ai(cos (pi^r^i +sin^t T ^.) 5 (4) 

i 

where is a four-component spinor at lattice site i 
and a is the lattice constant. The terms proportional 
to t and [i on the right-hand side of Q implement the 
term proportional to Sk in Eq. 0, while the terms on 
the second line accommodate the Zeeman splitting and 
the Rashba coupling. The terms on the third line cor- 
responds to the superconducting order parameter which 
is spatially varying. The normal and Rashba hopping 
terms couple only nearest-neighbour sites. We have al- 
lowed a spatially varying Rashba coupling and included 
a random potential Vi to study the effects of disorder. 
In the following we are working in the units where en- 
ergies are expressed in the units of the Rashba scale 
Er = ma 2 /2h 2 and the lengths are given in the units 
of the spin-orbit length l so = 2h 2 /ma (here a corre- 
sponds to the absolute value of the bulk Rashba cou- 



pling away from the junction). Typical values for these 
parameters InSb wires are Er ~ 50 jieV and l so ~ 200 
nm [lj. The disorder potential is modelled as a (dis- 
crete) Gaussian white noise with a correlation function 

(ViVj) = ^Y^-Sij characterized by the elastic mean free 
path l e and the Fermi velocity vp = ErI so /H. We are 
interested in the energy spectrum of a superconducting- 
normal wire-superconducting (SNS) junction where the 
left S segment is in the v = — 1 phase and the right one 
is in the v = 1 phase. Other relevant length scales in 
the problem are the superconducting coherence length 
£ = hvp/ A, the length of the normal segment Ln and 
the length scale ^r associated with the sign change of the 
Rashba coupling. In the following we model the variable 
Rashba coupling by aj = atanh [(j — N/2)a/^R] + a' 
with | a | > \a'\. Since we are working with a finite system 
of N lattice sites, the energy level crossings of topolog- 
ical origin in an infinite system are rendered to weakly 
avoided crossings that approach degeneracy as the effec- 
tive system size is increased. Motivated by experiments 
[lj, we consider wires with effective lenghts L > 10 l so (~ 
2-4 /im). 

Spectra of wires in the clean limit Vi = in the case 
B = (0, 0, B) are shown in Fig. |5J At cp = 0, 2tt chi- 
ral symmetry is unbroken and the Andreev states cor- 
responding to the MBS at the junction are degenerate 
within the resolution of the plot. The spectrum of a 
short junction Ln <C £, depicted in Fig. 2 a), has two 
bound states localized at the junction and qualitatively 
corresponds to a 7r-shifted spectrum of usual Majorana 
weak links [5] [131 E] • The inset illustrates how the spec- 
trum of MBS at the junction depends on the width of 
the Rashba domain wall £r. Longer junctions Ln > £ 
may have a number of states localized on the junction, 
as shown in Fig. [2]b), but also in that case spectrum is 
qualitatively shifted by tt. Qualitative features of spectra 
are insensitive to the relative lengths of Ln and £r. 

When magnetic field has a finite overlap with the spin- 
orbit field B x 7^ 0, chiral symmetry is violated and the 
strict topological distinction between v = ±1 phases dis- 
appears. This also implies that the (near) degeneracy 
of the MBS at (p = is lifted and moves to cpo ^ 0. 
Remarkably, as numerical calculations indicate, the vi- 
olation of chiral symmetry by finite B x is inefficient in 
lifting the degeneracy when the Rashba domain wall is 
not strongly asymmetric a' <C a. As shown in Fig. 3 
a), the offset is small (p>o ~ even for large symmetry- 
breaking fields B x ~ Er and significant asymmetry 
(a + a') /(a — a') = 1.5. Thus, the operation of the 
7r-j unction does not require a special fine tuning of pa- 
rameters. This behaviour is markedly different from the 
case of Majorana weak links where the sign of the Rashba 
coupling is constant. There the 5^-induced shift of the 
crossing point, taking place at odd multiples of tt when 
B x = 0, is large (cpo — 7r)mod27r ~ 1 for corresponding 
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Figure 2: a): Low- lying spectrum of a wire with a short 
junction (Ljy <C £) as a function of the phase difference, 
corresponding to parameters 5 = 4, A = 3, /i = and 
Z e = oo. The effective length of the wire is L = 10Z SO , 
£ R = 0.5 Z so and TV = 200. The flat mid-gap level is 
doubly degenerate, corresponding to the Majorana end 
states. The arcs correspond to the MBS at the junction, 
being nearly degenerate at ip = 0, 2tt. Inset: MBS 
spectrum at the junction for ^ = 0.1 Z so , l so and 2 Z so 
(outside to inside), b): Same as a) in the case of a long 
junction Ln = 4/ so = 12£. 




a) 7T 2vr b) 7T 2tt 

Figure 3: a): Spectrum of a wire with a short junction 
and asymmetric Rashba domain wall a' /a = 0.2 in the 
presence of large chiral symmetry breaking field B x = 1 , 

other parameters £? = 5, A = 3, /i = and l e = oo, 
L = 13Z so , £ R = 0.5 Z so and TV = 200. b): MBS spectrum 
of a short junction in the presence of disorder. The 
different curves correspond to different disorder 
configurations with l e = 2/ so , other parameters 5 = 5, 
A = 3, fi = 1, L = 10 / so , £r 0.5 Z so and N 200. 



values of B x . 

In Fig. 3 b) we have plotted 7r-j unction spectra in the 
presence of disorder. For a realistic disorder strength 
l e = 2 Z so the qualitative picture of clean wires survive 
although the spectra exhibit significant sample to sample 
fluctuations. Potential disorder affects the amplitude of 
the dispersion but cannot move the (avoided) crossing 
points. In the strong disorder regime l e <C l so the well- 
defined MBS localized at the junction are absent and the 
picture discussed in the clean limit breaks down. 

Discussion- In 2d semiconductor structures a finite 
Rashba coupling can be traced to structural asymmetry, 
resulting to a finite average electric field perpendicular 
to the sample plane. Moving electrons thus experience 
a magnetic field which couples to their spin. Structural 
asymmetry can be affected by asymmetric doping and ap- 



plying a gate voltage. A gate-induced sign reversal has 
been demonstrated in 2d systems [29J which is encour- 
aging for the present proposal. Semiconductor nanowire 
experiments show that the direction of the effective spin- 
orbit field is perpendicular to the wire and the substrate 
plain [TJ [30] . The sign of the Rashba coupling could be 
tuned, in principle, by local gating in wires also. Since 
the location of chemical potential needs to be controlled 
to achieve a topological segment, two independent gates 
are required. Although the double-gate tuning of the 
Rashba coupling in InAs wires has been experimentally 
realized [31] , combining this method with proximity su- 
perconductivity is a challenging task. 

Recently it was proposed that the Rashba coupling in 
nanowires could be emulated by spatially varying mag- 
netization [32-34J. As discussed in Ref. [34J, an effective 
Rashba coupling could be realized by a realistic array 
of nanomagnets in the vicinity of the wire. The effective 
Rashba coupling is proportional to the magnetic field gra- 
dient along the wire and may take values comparable to 
couplings in intrinsic Rashba wires [34J. In a configura- 
tion where the field is restricted to the sample plane, the 
sign of the Rashba coupling is determined by the rotation 
direction of the field along the wire [34J. If the magnets 
are stacked as shown in the inset of Fig. [l]b), the ro- 
tation of the magnetic field along the wire changes from 
clockwise to counterclockwise in the middle, effectively 
implementing the sign change of the Rashba coupling. 
Synthetic Rashba couplings are typically spatially mod- 
ulated which, according our numerical calculations, does 
not affect the operation of the tt-j unction significantly. 
Thus, engineering a topological 7r-j unction by employ- 
ing the method proposed in Ref. [34] should not require 
higher sophistication than realizing a wire with a cou- 
pling with a uniform sign. 

Conclusion- By controlling the sign of the Rashba cou- 
pling constant in superconducting nanowires, it is possi- 
ble to realize three topological phases protected by chi- 
ral symmetry. The interface between the two distinct 
topologically nontrivial phases supports two degenerate 
Majorana bound states and behaves as a topological tt- 
j unction when a superconducting phase gradient is ap- 
plied over the interface. This junction exhibits maximum 
supercurrent in the vicinity of vanishing phase difference 
and is robust against disorder and insensitive to moder- 
ate chiral symmetry breaking magnetic fields. 
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